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Abstract. In addition to revealing spherical harmonic degrees, I, of excited modes, pulsational 
amplitudes and phases from multicolour photometry and radial velocity data yield a valuable 
constraints on stellar atmospheric parameters and on subphotospheric convection. Multiperiodic 
pulsators are of particular interest because each mode yields independent constraints. We present 
an analysis of data on twelve modes observed in FG Vir star. 
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1. Introduction 

In a recent paper Daszynska-Daszkiewicz, Dzienibowski, Pamyatnykh (2003, Paper I) 
showed that photometric amplitudes and phases of pulsating stars are useful not only 
to identify the spherical harmonic degree, ^, but also for constraining models on stellar 
convection and stellar atmospheric parameters. 

To calculate theoretical values of photometric amplitudes and phases we have to make 
use of the complex nonadiabatic parameter, /, which gives the ratio of the local flux 
perturbation to the radial displacement at the photosphere. The / parameter is obtained 
from linear nonadiabatic calculation of oscillation in relevant stellar models. In the case 
of 5 Scuti star models, / is very sensitive to the treatment of subphotospheric convection. 
We also need models of stellar atmospheres for evaluation of the flux derivatives with 
respect to effective temperature and gravity and the limb-darkening coefficients. The 
most popular are Kurucz models (1998) but now there are alternative models available, 
e.g. Phoenix (Hauschildt et al. 1997) or NEMO. 2003 (Nendwich et al. 2004). 

In Paper I we proposed the method of extracting simultaneously i and / from multi- 
colour data and applied it to several 5 Scuti stars. In all cases the i identiflcation were 
unique and the inferred values of / were sufficiently accurate to yield a useful constraints 
on convection. In the follow-up paper (Daszynska-Daszkiewicz et al. 2004), we included 
radial velocity measurements, which improved significantly the determination of i and /. 
The very important feature of this method is that we can determine the i values in 5 Scuti 
stars avoiding major theoretical uncertainty concerning subphotospheric convection. 

Multimode pulsators are of special interest because each mode gives us independent 
constraints on stellar parameters. FG Vir is the most multimodal 5 Scuti pulsator. Af- 
ter last photometric and spectroscopic campaigns the total number of detected modes 
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Figure 1. Changes in photometric amphtudes and phases between 2002 and 2003. 



increased up to 48 (Brcger ct al. 2004, Zima et al. 2004). For twelve frequencies we have 
both two-colour Stromgren photometry {v,y), as well as radial velocity data. We apply 
our method to these modes using various models of stellar atmospheres. 



2. Observations 

Two recent photometric campaigns on FG Vir took place in 2002 and 2003. Spec- 
troscopic observations were carried out in 2002. In Fig.l we show how amplitudes and 
phases change from season 2002 to season 2003. We can see that for some of the modes, 
the differences are significant. This problem is more evident if we combine photometry 
and spectroscopy. In Fig. 2 we present diagnostic diagrams for the spherical harmonic 
degree, i, constructed from photometric and spectroscopic observables. The conclusion 
is that only contemporaneous observations can be used. 



3. Method of inferring i and / from observations 

The method is described in detail in Paper I. Here we give only a brief outline. 
The complex photometric amplitudes for a number of passbands, A, are written in the 
form of the linear observational equations 

V^ief)+£^e^A\ (3.1) 

where 

e^eYp{i,0), (3.2) 

, _ 1 ,mog(^fA|6^ 
""^-r^ aiogTeff (^-^^ 
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Figure 2. Diagnostic diagrams for the spherical harmonic degree, £, using combined 
photometric and spectroscopic data. 
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Having spectroscopic observations we can supplement the above set with an expression 
for the radial velocity (the first moment, A4i) 

GM 



ILjR 



(3.5) 



Symbols in Eqs. (3.3-5) have the following meaning, e is a complex parameter fixing mode 



amplitude and phase, i is inclination angle, and 



are limb-darkening-weighted 



disc averaging factors. J-{Tcs, g) is the monochromatic flux from static atmosphere mod- 
els. Mean atmospheric parameters, (Teg, \ogg^ [m/H]), enter through P^, £^ and through 
the disc averaging factors, which contain the limb-darkening coefficients. 

In the calculations reported in this paper, we used J-\{Tcs^\ogg,[m/ H\) from the 
Kurucz, Phoenix, and NEMO. 2003 models. The limb-darkening coefficients were taken 
from Claret (2000, 2003) for the Kurucz and Phoenix models, and from Barban et al. 
(2003) for the NEMO. 2003 models. 

Each passband. A, yields r.h.s. of equations (3.1). Measurements of radial velocity yield 
r.h.s. of equation (3.5). With data from two passbands we have three complex linear 
equations for two complex unknowns: e and (e/) (note that |e/| is the relative amplitude 
of the luminosity variations). The equations are solved by LS method for specified (. 
values. The £ determination is based on the behavior of minima of x^(^)- The inferred / 
values are of interest as constraints on stellar convection models (Paper I). 



4. Application to FG Vir 

EG Vir is a well known 5 Scuti variable located in the middle of its Main Sequence 
evolution. The basic stellar parameters, as derived from mean photometric indices and 
Hipparcos parallax, are: logTcg = 3.875 ± 0.009, logg = 4.0 ± 0.1. The mass estimated 
from evolutionary tracks is M = 1.85 ± O.lAf0. As was shown by Mittermayer & Weiss 
(2003), FG Vir has the solar chemical composition, thus in our calculations we adopt the 
standard metal abundance Z = 0.02. 
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Figure 3. The values of from the method described in the text as a function of £. 

Table 1. Possible identification of £ (x^ ^ 1-6) within the accepted T^s range 

v [c/d] our identification Breger et al. Viscum et al. Breger et al. 

from phot.&spec. (1995) (1998) (1999) 
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4.1. Identification of spherical harmonic degree, £ 

We applied the method to twelve frequencies of FG Vir. All of them have both photomet- 
ric and spectroscopic data. We use only observations from 2002 because radial velocity 
measurements are only from 2002. In the present application the radial velocity data 
are essential because we have data for only two photometric passbands and three is the 
minimum if we want to rely on the pure photometric version of our method. 

In Fig. 3 we plot as a function of £. The results were obtained using Kurucz models. 
We see that the discrimination of the £'s is sometimes better, sometimes worse. A unique £ 
identification is possible only for the highest amplitude modes. In Table 1 we compare our 
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Figure 4. Comparison of empirical values of nonadiabatic parameter, /, (dots with error bars) 
with the theoretical ones calculated for two values of MLT parameter, a. 



i identifications with earlier ones. Clearly, we are less optimistic than our predecessors. 
We rejected only the £ values leading to ^ 1-6- With this criterion we could assign a 
unique £ values only to three highest peeks. It is significant that £ > 3 is excluded in all 
twelve cases at the safe confidence level. 



4.2. Constraints on stellar convection 

In Fig. 4 we compare empirical values of parameters / deduced from the data with the 
values obtained from linear nonadiabatic calculation. Here we rely on a rather naive 
treatment of convection: mixing-length theory (MLT) and convective flux freezing ap- 
proximation. 

The empirical /'s depend only weakly on adopted values of TcS and L. The plotted 
values were obtained for M = 1.85Mq, logToff = 3.873 and logL — 1.156. These param- 
eters are consistent with mean photometric data for FG Vir and evolutionary models, 
and they lead for the three highest amplitude modes to the lowest obtain with our 
method. The empirical /'s are relatively insensitive also to the choice of £ which in some 
cases is ambiguous. For the plot we use £'s corresponding to minimum. The calculated 
/'s are even less ^-dependent. 

We can see not a bad agreement with theoretical values of the / parameter for models 
calculated with a — 0.0. This is good because our approximations are relevant in the 
limit of totally inefficient convection. 



4.3. Consequences of using different models of stellar atmospheres 

Our determination of £ and / requires atmospheric models. All our results so far were 
obtained with the use of Kurucz tabular data. In our comparison we use the data from 
Phoenix and NEMO. 2003 models. Such a comparison is important for assessing reliability 
of the quoted £ and / values. Models of stellar atmospheres are needed for evaluation 
coefficients I?^ , £^ and b^. The most important quantity derived from the models is the 
temperature derivative of the monochromatic fiux 
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Figure 5. as a function of effective temperature a tfie dominant mode, ui. Dots witli error 
bars sfiow the values derived from mean colors. 

which enters . This quantity depends most strongly on T^s and the same must be 
true for the minimum x^- This suggest that the best value of Tcff in the same as the best 
value of £. 

We found that relying on different atmospheric models, the £ identification based on 
minima is unchanged but the minimum values of in the accepted T^g range were 
much larger than the ones obtained with the use of Kurucz data. We then relaxed the 
temperature constraint and considered models in a wide range of Toff. We use the fit of 
z/2 to the radial fundamental mode and evolutionary tracks to derive corresponding value 
of log f;. We encountered two problems, which are illustrated in Fig. 5: 

• acceptable minima of for Phoenix and NEMO. 2003 models were found outside 
the acceptable Tcff range 

• in Kurucz and Phoenix models the temperature flux derivatives, ax, are non-smooth 
which gives bad shape of x^iTes) 



5. Conclusions 

A unique identification of spherical harmonic degree, £, of excited modes without a 
priori knowledge of complex parameter /, which links the surface flux variation to dis- 
placement, is possible. However, in the case of low amplitude modes, more accurate 
observations and measurements in more passbands are needed. 
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The values of / inferred from data for FG Vir are crudely consistent with models 
calculated assuming inefficient convection {a « 0.0). We showed that there is a prospect 
for constraining Toff from the pulsation data. However, for this application wc need more 
accurate derivatives of the fluxes in various passbands with respect to T^g than the 
derivatives we may calculated from available tables of atmospheric models. 

There are indeed various applications of amplitude and phase data from photometric 
and spectroscopic observations of pulsating stars. Multimode pulsators are advantageous, 
but the mode amplitudes have to be determined with a high accuracy. Data from pho- 
tometric observations in more than two passbands are very much desired. 
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